function [Polos, Zeros] = correction(porcentaje, WN, Polos, Zeros, KK, aa)

            sys = KK*tf(poly(Zeros),poly(Polos));

            w=logspace(log10(1),log10(WN*1.1),10000);
            [MAG]=bode(sys,w);
            mag=-squeeze(20*log10(MAG));
            encontrado=(find((mag>aa-0.001),1,'first'));
            if encontrado>0
                Wacalc=w(encontrado);
                correccion=(WN-Wacalc)*porcentaje/Wacalc;
                
                Zeros=Zeros*(1+correccion);
                Polos=Polos*(1+correccion);                
            end


end